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Binary interactions with high accretion rates onto main sequence stars 

Sagiv Shiber, Ron Schreier and Noam Soker^ 

ABSTRACT 

Energetic outflows from main sequence stars accreting mass at very high rates 
might account for the powering of some eruptive objects, such as merging main se¬ 
quence stars, major eruptions of luminous blue variables, e.g., the Great Eruption of 
Eta Carinae, and other intermediate luminosity optical transients (lEOTs; Red Novae; 

Red Transients). These powerful outflows could potentially also supply the extra en¬ 
ergy required in the common envelope process and in the grazing envelope evolution of 
binary systems. We propose that a massive outflow/jets mediated by magnetic fields 
might remove energy and angular momentum from the accretion disk to allow such 
high accretion rate flows. By examining the possible activity of the magnetic fields 
of accretion disks we conclude that indeed main sequence stars might accrete mass at 
very high rates, up to ~ yr“^ for solar type stars, and up to ~ IMq yr“^ for 

very massive stars. We speculate that magnetic fields amplified in such extreme con¬ 
ditions might lead to the formation of massive bipolar outflows that can remove most 
of the disk’s energy and angular momentum. It is this energy and angular momentum 
removal that allows the very high mass accretion rate on to main sequence stars. 

Keywords: Stars: stars: AGB and post-AGB - stars: jets - (stars:) binaries (including 
multiple): close 


1. INTRODUCTION 


Evidences are acquired in recent years on the need to allow compact objects to accrete mass 
at rates exceeding the Eddington limit by about one to two orders of magnitude. Extreme cases of 
super-Eddington accretion rates is the magnetized neutron sta rM82X-ll that accr etes from a mas¬ 
sive companion at a rate of ~ 100 times the Eddington limit (IBachetti et al.ll2014h. an d two active 


galactic nuclei (AGNs) with 200 — 300 times the Eddington limit (IDu et al. 


2010). Other ex- 


amples include accretion onto a main sequence star during mergers of stars ( Soker & Kashil 2011 ) 


during periastron passages of main sequence stars near asymptotic giant branch stars (IStaff et al. 
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2016r) . an d outbursts of some ma ss ive stars that migh t be power by accretion onto a main sequence 
star (e.g., iKashi & SokenbOlOallbl : iKashi et alJl2013b : a general name for these type of objects is 
intermediate luminosity optical transients (ILOTs; also termed Red Novae or Optical Transients). 

In many cases, super-Eddington or Sub-Eddin gton outflow kinetic p ower is larger than th e radia¬ 
tion luminosity, e.g., the AGN HS 0810-1-2554 ( Chartas et al. 2014 ). Chartas et al. ( 2014 ) suggest 
that this implies that magnetic activity is likely behind the acceleration of the outflow. 

Despite the observational hints for high mass accretion rates onto main sequence stars, it is 
comr nonly assumed that main s equence stars cannot accrete mass at high rates, above ~ IO^^Mq yr“^ 
(e.g.. lHjellming & Taamlll99ll) . Our main theme in the present study is to introduce the idea that 
main sequence stars might accrete mass at very high rates. This has implications to transient 
events, to the common envelope evolution, and to the newly proposed grazing envelope evolution. 
We approach this challenge by building on a strong magnetic activity in the accretion disk, in¬ 
cluding an efficient dynamo and magnetic field reconnection. We deal only with main sequence 
stars, and with accretion rates of up to ~ O.OIMq yr“^ for low mass main sequence stars, and up to 
~ IMct, vr~^ for yery hig h mass main sequence stars. In the merger model for the lEOT V838 Mon 
dSoker & Tvlendal l2003h . for example, a companion of O.SM© was destroyed on a star of 


~ GMq. It is assumed that several percents of the companion mass was accreted within several 
months; the r est inflated a l a rge en velope. This leads to an accretion rate of ~ 0.01 — 0.1 Mq yr“^. 
We note that iNandez et al.l (l2014h simulate a merger process with parameters appropriate for the 
lEOT VI309 Sco, and find the formation of an accretion disk in the synchronized case. They did 
not find jets, as expected if magnetic fields are not included. As well, in the non-synchronized case 
the merger proceeds without the formation of an accretion disk. Namely, it is possible, according 
to their results, that merger and high accretion rate does not necessary requires an accretion disk. 

It is commonly agreed upon that jets are produced at the center o f accretion disks and that 


magn etic fields are a crucial ingredient in their launching process (e.g., lEivioll201ll: IPudritz et al. 


20121). Many models for the formation of astrophysical massive outflows are b ased on the oper¬ 
ation of l arge scale magnetic fields within and around the accretion disk (e.g., IZanni & Eerreira 
20131 and iNarayan et al.ll2014l and references therein). By large-scale we refer to magnetic fields 
with coherence scale larger than the radius o f the disk at the consid ered location. In the ”X-wind 
mechanism” introduced by IShu et al.l (Il988b and IShu et al.l (Il99ll) the jets are launched from a 
narrow region in the magnetopause of the stellar field. A different group of models also use open 
radial magnetic field lines and are based o n an outflow froni an extended disk region, but they dp 
not rely on the stellar magnetic field (e.g., Konigl & Pudritz 2000 : Shu et al. 2000; Eerreira 2002 : 
Krasnopolskv et al. 2003; Eerreira & Cassa 2004 ). 


There are several arguments that point to problems with models for launching jets that are 
based only on large scale ordered magnetic fields. (1) Precessing jets. In several young stellar 
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objects (YSOs) the jets precess on a time scale of 100 years (e.g. lYbarra et alJl2006l) . A large 
scale magnetic field cannot change its symmetry axis on such a short time. (2) Collimated jets 
in planetary nebulae (P Ne). There is a highly collimated clumpy double-jet in the PN Hen 2-90 
dSahai & Nvmanll2000h very similar in properties to jets from YSOs that form Herbig-Haro (HH) 
objects (dense blobs along jets launched by YSOs). The source of the accreted mass in PNe is 
thought to be a companion star. In such a system large scale magnetic fields are not expected. (3) 
No jets in DQ Her (intermediate polars) systems. DQ Her systems are cataclysmic variables where 
the magnetic field of the accreting white dwarf (WD) is thought to truncate the accretion disk in its 
inner boundary. This magnetic field geometry is the basis for some jet-launching models in YSOs 
(e.g., IShu et al.lll99lh . However, no jets are observed in intermediate polars. (4) Some of these 
models don’t take into account the stro ng dynamo that is expected to operate in accretion disks. 
In few cases, e.g., iHujeirat et al.l (l2003t) . the modification of the large scale magnetic field in the 
disk is considered. Such dynamos substantially modify the structure of the rnagneti c field in the 
accretion disk and lead to a behavior closer to that in active stars. (5) Ferreiral (12013h summarizes 
accretion models of classical T Tauri (TTS) stars, and concludes that steady state models, such as 
X-winds or Accretion Powered Stellar Winds, cannot spin down a protostar. The balance between 
the accretion torque (working so as to spin up the star) and the magnetic torque (working so as 
to slow the stars’ rotation) should lead to a spin up on a time scale smaller than the disk lifetime, 
whereas observations show that the rotation period of TTSs is almost constant, at roughly 10% 
of their break-up speed. Instead, unsteady models, such as Magnetospheric Ejection are the best 
candidates for removal of extra energy from the disk. The model suggested in this study is a type 
of unsteady model. 

In the present paper we take the view that the yery high accretion rat es onto main sequence 
stars in ILOTs, in the range of ~ 10““^ — IMq yr“^ (IKashi & Sokenl20I0bh . require the formation 
of thick accretion disks. The accretion rates in our model are high, and the disks are thick enough 
so as to sustain turbulent eddies for long durations. Those eddies, combined with the differential 
rotation of the disk’s plasma, generate the dynamo magnetic fields. Strong fields will eject mass 
and momentum out of the disk without the need of inward diffusion. In section |2] we build such 
accretion disks that we will use in section [3] where we study their magnetic field properties and 
activity. We will take the magnetic dynamo to be similar in many respects to that of a stellar 
dynamo. The merit of comp arison to the dynamo operating in the Sun is demonstrated in a recent 
paper by lBugli et al.l (120141) . They study dynamo inside a thick disk around a Kerr black hole and 
find the evolution of the magnetic field to qualitatively occur in the same fashion as in the Sun. 
Our summary is in section HI 
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2. BUILDING THE DISK 
2.1. A thick disk 


We here construct an accretion disk with a high accretion rate w here part of the energy is 
carried away by jets. High accretion rate disks were studied before, e.g.. lAbramowicz et al.l (Il980t) 
for an accreting BH, as we ll as accretion disks wh ere heat is removed to h eat the corona (e.g., 
Svensson & Zdziarskilll994t) . or to an outflow te.g.. iKuncic & Bicknellll2006t) . To present a com¬ 
prehensive picture we repeat the construction of the disk presented in these studies, and present it 
in a way that will serve our goals. 


disks (Shakura & Sunvaev 1973 


Svensson & Zdziarski (Il994h modified the st andard Shakura-Sunyaev model for accretion 


Frank et al.ll2002h by allowing a major fraction, /, of the power 


released from the accreted disk matter to be transported and dissipated in the corona. We replace 
their (1 — /) factor by the parameter e that is the fraction of the energy liberated by the accretion 
process that goes to radiation 


AaT^ _ 3GMM 
3r SttR^ 



( 1 ) 


In our study the rest of the energy is channelled to the outflowing gas, while the contributions to 
the pressure from this channels, e.g., magnetic fields, are neglected. The other symbols have their 
usual meaning: M is the mass of the star, M is the mass accretion rate, R is the radial coordinate in 
the disk plane, i?* is the stellar radius, r is the optical depth, and a is Stefan-Boltzmann constant. 
We study accretion disks with very high mass accretion rates of ~ 10“'^ — IMq yr“^ in which the 
primary star is a main-sequence star with a radius of i?* ~ 10^^ cm. This can cause the opacity 
to be dominated by electron scattering processes rather than Kramers opacity induced by free-free 
emission, and also the pressure to be dominated by radiation pressure rather than gas pressure. 

In the gas pressure-dominated regime we take only gas pressure, starting with the Kramers 
opacity. We will use a molecular weight of /r = 0.615 throughout the paper. When the parameter e 
is included, we find the following expressions for the radial dependance of the disk height H, the 
central temperature Tc, and the radial velocity of the disk vr, respectively 

H/R = (/^M24) , (2) 

T, = 2 X K, (3) 

and 

VR = -3.5 X V'/' cms-\ (4) 
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where a is the viscosity disk parameter which assumed to be radius-independent, and we defined 
the dimensionless variables 

M 

m = Ri2 = 


M24 = 


Me 

M 



g \i !-> ( ^ 

For these typical values the disk is no longer thin, and the radial velocity has the same order 
of magnitude as the Keplerian velocity. In addition, the high temperature suggests that electron 
scattering processes and radiation pressure might be important. 

Taking the pressure to be radiation pressure and taking electron scattering opacity give the 
following disk properties 

H/R = 1.656 


'"M24 

)R 

- 1/4 

m 


-1 
12 5 


Rl 


12 


1/8 


K, 


( 6 ) 

(7) 


and 


vr =—A.Q5 X 10'^ cms (8) 

The properties are more sensitive to the parameter e, which can str ongly affect the thickne ss and 
the radial velocity of the disk. As noted in previous works (e.g., ISocrates & Davisll2006t) . it is 
possible to build thin disks even for very high accretion rates by taking 6 -C 1. Another important 
difference is that the disk height over radius is a decreasing function with radius where in the gas 
pressure dominated case it is an increasing function with radius. 

In any case, in order not to required extremely small value of e, and demanding that H is not 
much larger than r, we find from the expressions for if/r that the accretion rate in our treatment is 
limited. To an order of magnitude this limitation on the accretion rate is ~ yr“^ for solar 

type stars, and up to ~ IMq yr“^ for stars of 30 and more solar masses. 


2.2. Relevant time scales 


There are several relevant time scales. The dynamical timescale (or the Keplerian time divided 
by 27r) is given by 

^ ^ (9) 

Vcj) 


This time scale is about equal to the sound crossing time along the disk height, which is about the 
dynamical time scale to reach equilibrium along the 2 ; direction (perpendicular to the disk plane) 

H 


t 


c. 


( 10 ) 
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where Cg is the sound speed. The viseous time seale is given by 

R 

_ 

Vise ~ ~ 1 

V Vr 


where v = aCgH is the turbulence viscosity coefficient. The thermal timescale is given by 

heat content per unit area 
dissipation rate per unit area D {R) 


( 11 ) 


( 12 ) 


Where D {R) is the dissipation due to differential viscous torque and S is the surface density of 
the disk. 


In the Radiation dominated regime it can be shown that 

, H tdis 

tth = £—T = e—, (13) 

c 3 

where tdis = ‘iHr/c is the photon diffusion time from the disk’s mid-plane out. From equation 
(fT3]) we can interpret £ as the probability to lose energy (heat) via photon diffusion (hence radiation) 
rather than outflow kinetic energy. Approximately, as the diffusion time becomes longer a larger 
fraction of the energy is transported by other than radiation channels. In our model these are 
magnetic fields and gas outflows. We note that when magnetic fields remove energy much faster 
than photons do, then the disk might be geometrically thin, but still the radiation carries only a 
small fraction of the energy. 

Let us dwell on this point. The viscosity in the disk and magnetic activity dissipate energy. 
But in our regime of very high accretion rate, the photon diffusion time is very long. Most of 
the dissipated energy that starts as thermal energy, rather than being emitted by the disk, will 
lead to internal motion, such as turbulence that amplifies magnetic fields, and to gas buoyancy 
outward. Near the surface, in particular due to magnetic field reconnection below the photosphere 
(see section |3]), the energy is channelled to winds or jets. The photons have no time to diffuse, and 
most of the energy is carried by the outflow, i.e., kinetic energy. The acceleration occurs in many 
sporadic impulsive events resulting from magnetic field reconnection. Most of the reconnection 
events are much shorter than the dynamical time, and hence the mass and energy they eject merge 
to one continues outflow. 


For the treatment above to hold, the inflow time must be longer than one rotation time of 
material in the disk 

Un > 4ep = 27rt0. (14) 

This condition can be cast to height over radius condition. With the timescales definitions fin = 
R/vr, = H/cs and using the relation vr = (3z//2i?) one gets 


1 3 

— ^ a— 
271 ^ 2 




( 15 ) 
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The area in the e — M plane where this inequality holds is termed the ‘safe-zone’. For others values 
of a, R^, and M the safe-zone changes. 


3. MAGNETIC ACTIVITY 


3.1. Dynamo 

The high accretion rate is associated with very large Reynolds flows and therefore the transi¬ 
tion to turbulence is inevitable. The differential rotation in the disk is responsible for the generation 
of the toroidal magnetic field, the fl-effect, and the combined action of helical turbulence and dif¬ 
ferential rotation generates the poloidal field, the a-effect, and so forth. This is the bas ic of the 


a-fl dy namo. We follow the treatment of the magnetic field amplification as conducted bv IPudritz 
Jl981ah . 


We study flow regimes with very high Reynolds numbers. iNath & Mukhopadhyayl (120151) 
have shown recently that for high Reynolds numbers, as is t he case for realistic astrophysical accre¬ 
tion d isks, nonlinearity is achieved by transient growth ( e.g., Umurhan & Regevl2004 ^ Umurhan et al, 
2007 ) rather by the Magnetorotational Instability (MRI). Nath & Mukhopadhvav ( 2015 ) also showed 
that when the transient growth modes get to play, the magnetic fields in crease much mor e than what 
the pure MRI model gives. For that, although the results obtained by IPudritz (Il981ah might not 
be applicable when MRI dom i nates, i n the regime studied here the transient growth dominates and 
we use the results bv IPudritz I (Il981al) . 


The dynamo action is described by the mean field theory, in which the magnetic field and the 
velocity are decomposed to mean and fluctuating components, 

Bt = B + b, Ut = U + u. (16) 

The induction equation for the fluctuating component of the magnetic field is written under the so 
called ‘first order smoothing theory’ as 

db 


dt 


= Vx [U y. b -\- u X B] + b 


(17) 


The last term comes from Ohmic dissipation of the fluctuating component of the current, given 
by i = (c/47r)V x b, and reduces the magnetic field. The first term on the right hand side might 
lead to convective amplification of the magnetic field, the ‘a-dynamo’. From the assumption that 
the Ohmic dissipation rate of fluctuating magnetic energy per unit volume eoh = 
equals the rate at which turbulenc e pump energy back Ctr = {tit/IDB'^/Ati, where lu is the largest 
turbulence scale, IPudritz (Il981al) derived the magnitude of the fluctuating magnetic field 


b^ = ^B^. 
r] 


( 18 ) 
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The magnetic diffusivity r] is given by the expression 

T] = c^/dvrcTs = 2.95 x 10^^ In A cm^/s, (19) 


where as is the Spitzer conductivity, T is the temperature, and In A is 
r]T is the turbulent diffusivity for the mean magnetic field B, 


of order 10 (see Huba 112013 ). 


= M^H^/tK. ( 20 ) 

Here H is disk height, tx the Keplerian orbital period, and Mt is the turbulent Mach number 
defined as 


Mt = u/cs, 

with u the r.m.s. turbulent speed, and Cg is the speed of sound. 

Substituting equation ^ into (l20l) . and dS]) into (fT9l) we find the ratio 


Mt\" / £ \ V4 


= r]T/v = 3.2 X 10 ^^ ( ^ ) ( ^ 


X 


a \ - 1/2 / m 
(U 


Rh 


0.5 J VO.l 

^ 4 • \ 3/4 

/"M 24 ) 


( 21 ) 


( 22 ) 


There are several large uncertainties in the derivation of equation (l22l) . Nonetheless, its im¬ 
plication is clear: the fluctuating magnetic field dominates. Any initial large scale, of the order of 
the disk radius and more, magnetic field becomes negligible in the disk. At these high accretion 
rates large scale magnetic fields are not expected to play a significant role in the ejection of jets or 
disk winds. We should look for the role played by the fluctuating magnetic field, that is amplified 
by the dynamo, in launching the jets and winds. 

The ratio of the fluctuating magnetic pressure to the total pressure is defined by 

A = b'^/A'irP, (23) 


so that the mean magnetic field is 


Taking P 


B = \/MTp~qJ^. 


pel. 


Pudritz & Fahlmanl (119821) estimated a very large value A of 


10^/2 < y < loA 


(24) 


(25) 


As we are discussing accretion into main sequence stars where we apply the disk structure neglect¬ 
ing magnetic pressure (section [2l), we limit the ratio of the magnetic to gas pressure to A < 1. As 
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well, we take a eonservative approaeh in this preliminary study. We derive the following value for 
the magnetic field 


\b\ =2 X 


£ 


X 


m 

74 


0.1 

7/16 


- 1/40 




24 


a 

(U 

17/40 


- 9/20 


G. 


(26) 


Two comments are in place here regarding the numerical value given in equation (1261) . (i) 
In solar spots the magnetic fields reach equipartition, between the magnetic and thermal pressure, 
for magnetic fields of ~ 2000 G. The magnetic fields in the accretion disk studies here are about 
one order of magnitude stronger. This implies that the magnetic activity in the accretion disk is 
extremely strong, with magnetic pressures about two orders of magnitude above the photospheric 
pressure of the accreting star. It should be noted that for the very high accretion rates used here 
> 10“^Mq yr“^ and the low value of e ~ 0.1, both the density and temperature inside the disk are 
much higher than the corresponding values on the stellar photosphere. The pressure inside the disk 
is therefore higher by more than two orders of magnitude than that on the stellar surface, and the 
strong magnetic field is still below the thermal pressure inside the disk, (ii) The scaling in equation 
(l26l) takes A = 1. This is applicable to the case where the magnetic fields reaches equipartition 
with the thermal pressure in the accretion disk, as is the situation in solar spots. Above the solar 
surface the magnetic pressure can be much larger than the thermal pressure. Similarly, it is quite 
possible that in the accretion disk we study here the magnetic pressure will be larger than the 
thermal pressure, and rather be equal to the ram pressure of the Keplerian motion in the disk pv'^. 
This is the case for A > 1, as taken hv IPudritz & Fahlmanl(ll982h . Such magnetic fields will make 
the proposed scenario more feasible even. 


3.2. Magnetic power 

We show now that most of the energy could potentially be carried off by magnetic activity 
rather than by radiation even when the magnetic energy in the disk is smaller than the thermal 
energy. 

We assume that most of the magnetic activity takes place just inward to a radius Rb, not too 
far from the stellar surface. We neglect further amplification of the magnetic field inward to this 
radius. The mass accretion rate is given by M AttRbH pv^. We define the ratio of the inward 
radial velocity to the Keplerian velocity 5 = vr/v^ < 1, and the kinetic energy density in the disk 
Cd — pv‘^/2. The accretion power of material crossing the radius Rb inward, i.e., the energy flux 
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through radius RbAs given by 

Li,cc = ~ ‘27r6HRBV^ed. (27) 

We take the time to amplify and release magnetic energy to be the dynamical time given by equa¬ 
tion (fTOl) tb ^ ^ H/cs ~ Rb/v^. With an energy density = b'^/Sir the magnetic power 

within radius Rb is Lb ^ 2TrR^HeB/rB- The ratio of magnetic energy density to kinetic energy 
density is then 

^ « sR. (28) 

-^acc 

We recall that in deriving equation (l28l) we assumed that there is a global steady state, and hence 
the amplification time of the magnetic field is equal to the time scale to remove magnetic energy by 
reconnection. Hence, for a given power, shorter amplification and reconnection time scales imply 
weaker magnetic fields. Alternatively, if the ratio of the magnetic energy density to the kinetic 
energy density is some fixed value, like unity, then the magnetic power increases with decreasing 
time scales. 

As 5 < 1, the above simple estimate shows that the magnetic energy can be lower than the 
kinetic energy in the disk, and yet the magnetic activity might carry a large fraction of the accretion 
energy. For that the magnetic energy should be vigorously amplified by a vigourous turbulence in 
the disk. We assume that the high accretion rate leads to this behavior of the disk. 


3.3. Reconnection and outflow 


In this subsection we speculate that the magnetic activity might lead to an energetic outflow. 

We define three general coherence scales of the magnetic fields. Small scale defined as much 
smaller than the disk scale hight Ib ^ H. Moderate scale for which the coherence scale is of the 
order, but smaller, than the disk scale hight, 0.3i7 < Ib < H. Finally, the large scale fields is 
with Ib ^ r > H. In the present study magnetic fields with moderate scale, about equals to the 
turbulence scale, dominate. We give no role to large scale magnetic fields (unlike many other jets’ 
models). 


As magnetic flux accumulated and the tension increases due to the dyn amo action, mqderate - 
scale turbulence trigger local explosive events of magnetic reconnection (ILazarian et al.ll2015h . 
Those events are sporadic and decay as the extra magnetic energy is releas ed. At any given radius 
R the largest turbulence scale, R, is given hy R = Mf H (IPudritz Ill981at) . and for high turbulent 
Mach numbers such as Mt ~ 0.5 expected for the presently studied accretion flows (section^, the 
vertical scale can be populated by only 2 such large eddies. The large magnetic flux tube will also 
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buoy quite fast, relative to small flux tubes, outward from the disk, allowing the replenishing time 
of the magnetic fields to be of the order of the dynamical time. 


There are two key processes resulting from the large flux tubes. These are magnetic torque 
and reconnection in scales not much smaller than the disk scale height. In accretion disk models 
where the large scale magnetic fields extended beyond the disk, there is a magnetic torque on the 
disk from gas outside the disk, mainly a disk wind. When there is no such a field, torque acts within 
the disk. When the flux tubes are small relative to the radius, the torque is local and play no large 
scale role. In our situation the flux tubes are within the disk, but they are very large. Therefore, 
we speculate that they might transfer angular momentum and energy from the internal parts of the 
disk to near the surface of the disk. Near the surface magnetic field are likely to reconnect and 
accelerates the wind/jets. This process further reduces the amount of energy that is radiated from 
the disk, hence lowering the value of e. 


The vigourous turbulence facilitates magnetic field reconnection (ILazarian et al.ll2015h . We 
propose therefore, that magnetic field reconnection will take place not only above the disk, as in 
the Sun, but also below the disk surface. The outcomes of the reconnection process of the large 
flux tubes might be as follows. 


(1) Mass outflow. Like in the Sun, once reconnection occurs between two adjacent magnetic 
flux tubes plasma flows along reconnected magnetic field lines, in two opposite directions. In our 
setting we take that a large fraction of the reconnection events takes place below the photosphere. 
Although initially a large fraction of the dissipated energy is thermal, before photons have time 
to diffuse out the gas is accelerated to form an outflow. Each reconnection event is like a small 
explosion. Much as in supernova explosions where photons diffusion time is very long and most 
of the final energy is kinetic with only a small fraction emitted as radiation. 


(2) Huge magnetic arcs. The reconnection process implies more than direct mass ejection, 
and migth lead to the build up of huge magnetic arcs above the disk. While the original flux tube 
emerging from the disk are expect to reside at about the same radial distance R, we speculate 
that the new arcs formed by reconnection can have their two foot-p oints at different radii, much 
as trans-equatorial magnetic loops in the Sun (ITadesse et al.ll2014t) . We term these speculative 
structures trans-radial loops. The trans-radial loops experience very large shear, and might lead to 
further amplifi cation of the magnetic ene rgy. This process will be studied in the next paper in the 
series (see also lHeyvaerts & Priestil 19891) . 


(3) Flares. Reconnection of such huge arcs might lead to extremely energetic flares. Such 
an energetic fl are was speculated to be the caus e of the 2009 Novem ber outburst of Aqu i la X-1 
( Soker 2010 ). de Gouveia Dal Pino et al. ( 2010 ). Khiali et al. ( 2015 ) and Kadowaki et al. ( 2015 ) 
discuss such rapid reconnection processes in relation to different properties of accretion disks of 
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microquasars. We aim here at aeeretion disks around main sequenee stars aeereting at a very high 
rate. We expeet the wind from the disk to obscure radiation from sueh flares. 


4. SUMMARY 


We deseribed a erude aeeretion flow setting to allow main sequenee stars to aeerete mass at 
very high rates. The key to allow sueh high aeeretion rates is that most of the energy liberated in 
the aeeretion proe ess, and the angular rnqment um earried by the inflowing gas, are removed by 
massive outflows. ISvensson & Zdziarskil (119941) used the removed energy to heat the eorona, and 


Kuneie & Bicknelll (120061) produeed jets by magnetie torque, both in aeeretion disks onto blaek 


holes. We here studied aeeretion onto main sequenee stars, and the removal of energy and angular 
momentum by jets that are formed from sporadie reeonnection events of magnetic fields that are 
amplified by a strong dynamo. While in previous papers of our group it was just assumed that main 
sequence stars might accrete at very high rate, in the present study we have outlined the seenario 
for that. The main result is that the removal of angular momentum and energy eannot be aehieved 
without magnetie fields that are effieiently amplified by the operation of a dynamo in the dense 
aeeretion disk. 

As was done in the past (e.g., Svensson & Zdziar^ll994 ) we deseribed a disk strueture where 
only a fraction £ -C 1 of the energy liberated by aeeretion is ehannelled to radiation (seotion[2l). 
From this crude structure we eould eonstrain the properties of the disk around main sequenee stars. 
For the parameters used here we found a thiek aeeretion disk dominated by thermal pressure (eqs. 
[2]- HI). To still be regarded as an accretion disk in the usual meaning, the inflow time should be 
longer than the orbital time. 


Using the disk properties, we speeulated on possible magnetie aetivities that might lead to 
rapid mass ejeetion from the disk. We presented the expeeted dynamo amplifieation of the mag¬ 
netie field, and eoneluded that a very strong fluetuating magnetie field is expeeted in the disk (eq. 
[22l) . We then argued ('section [3^ that it is possible for the magnetie power to be very large without 
altering the general disk structure. We also note that for geometrieally thin disks, for which 5 1, 

the removal of energy by magnetic fields ean be effieient even for <C e^. Sueh a thin disk might 
be obtained if the disk is radiation dominated, as seen from equation (© for e 1. This outlined 
seenario for allowing high mass aeeretion rates onto main sequenee stars, although speeulative in 
some parts, is our main new result. We hope it will eneourage relevant models of astrophysieal sys¬ 
tems, sueh as the eommon envelope evolution, to eonsider main sequenee stars that aeerete mass 
at very high rates, as well as motivate numerieal studies of the aeeretion flow. Sueh 3-dimensional 
numerieal studies should inelude magnetie fields and be of high spatial resolution to resolve the 
turbulence in the aeeretion disk. 
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Some key ingredients of the proposed scenario should be emphasized. (1) In the proposed 
model magnetic fields with moderate coherence scales dominate. By moderate scale we refer 
to scale of 0.3H < Ib < H, where H is the disk scale-height. Large-scale magnetic fields of 
Ib ^ H play no role in our model (unlike in many other jets’ models). (2) Although magnetic 
energy is dissipated locally to heat by reconnection, in our setting most of this energy ends up 
as kinetic energy of the outflow. The reason is that because of the very high mass accretion and 
outflow rates, the flow is optically thick. Therefore, before photons have time to diffuse out the gas 
expands. Most of the thermal energy goes to adiabatic expansion rather than to radiation (much as 
in supemovae explosions). 


The outcomes of strong magnetic field reconnection below the photosphere, discussed in sec¬ 
tion [33l br ought us to propose the e xistence of trans-radial magnetic field lines, similar to those 
studied by (IHevvaerts & Priestill989h . Namely, huge magnetic arcs above the disk that have their 
two foot-points at different radii. These trans-radial magnetic fields might lead to further magnetic 
held amplification. We expect trans-radial loops to play a major role in the magnetic activity of 
the boundary layer, where the disk angular velocity substantially decreases to match the angular 
velocity of the star. The dynamics of trans-radial loops is the subject of a forthcoming paper. 


Our main conclusion is that main sequence stars might accrete mass at very high rates, up to 
IO^^Mq yr“^ for solar type stars, and up to ~ IMq yr“^ for very massive stars. Such accretion 
rates could potentialy account for some outbursting objects that have their luminosity in the bulk 
region between novae and supernovae and t hat are thought to be p o wered by high accr etion rates 
onto main sequence stars in binary systems (IKashi & Sokerll2010blJbl:lKashi et al.ll2013h . As well, 
such high a ccretion rates mig ht take place in main sequence companions in eccentric orbits around 
AGB stars dStaff et al.ll2016h . A general name for these types of objects is intermediate luminosity 
optical transients (ILOTs; also termed Red Nova e or Optical Transients) . Examples of specific 
objects include stellar merger, s uch as V838 Mon dSoker & Tvlendall2003h . the nin eteenth century 


Great Eruption of Eta Carinae dKashi & Sokerll2010ali rDrogenitors of some PNs dSoker & Kashi 


2012h. such as KiPn 8 dBoumis & Meaburnll2013f) . and the pre-explosion outbursts of SN 2009ip 


dSoker & Kashill2013l) . 


High accretion rates onto, hence high outflow rates from, main sequence stars might supply 
extra energy to remove a common envelope where the secondary is a main sequence star. When the 
envelope removal by jets is sufficiently efficient as to remove the entire envelope residing outside 
the secondary star, no real co mmon envel ope phase is established. Instead, a grazing envelope 
evolution (GEE) commences dSokeill2015h . The GEE might nicely explain the observation that 
the kinetic ene rgy of the bipolar out flow from the binary system HD 101584 exceeds the released 
orbital energy dOlofsson et al.ll2015h . 


We can summarize our study by stating that the possibility for main sequence stars to eject 
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winds/jets at extremely high rates, as we argued here, opens a rieh variety of proeesses that eould 
potentialy aeeount for some properties and some evolutionary puzzles of binary stellar systems. 

We thank an anonymous referee for very detail report and many valuable suggestions. 
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